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Abstract. We report an analytical investigation of the force and torque acting upon a 
magnetic dipole placed in the vicinity of a moving electrically conducting nonmagnetic 
plate. This problem is relevant to contactless electromagnetic flow measurement in 
metallurgy and extends previous theoretical works (Thess et al. Phys. Rev. Lett., 
96(2006), 164501; New J. Phys. 9(2007), 299) to the case where the orientation of 
the magnetic dipole relative to the plate is arbitrary. It is demonstrated that for 
the case of low magnetic Reynolds number the three-dimensional distributions of the 
induced electric potential, of the eddy currents and of the induced magnetic field can 
be rigorously derived. It is also shown that all components of the force and torque 
can be computed without any further approximation. The results of the present work 
serve as a benchmark problem that can be used to verify numerical simulations of more 
complex magnetic field distributions. 
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1. Introduction 

When an electrically conductive macroscopic body moves in the nonuniform magnetic 
field B created by an external source, for instance a permanent magnet, then eddy 
currents j are induced in the body. These eddy currents create the Lorentz force 
F = j X B. The force is directed in the direction opposite to the direction of the body's 
translation. Moreover, the eddy currents induce an additional magnetic field b which 
interacts with the permanent in such a way as to create a force and a torque acting 
upon the magnet. Thus, by measuring these forces and torques one can determine the 
parameters of the movement such as the velocity or the direction. This principle is 
embodied in a contactless electromagnetic fiow measurement technique called Lorentz 
force velocimetry [8], [10], [9], [1], [2], [6], which permits fiow measurement in hot and 
aggressive fiuids like liquid aluminium or molten steel. When developing measurement 
systems embodying Lorentz force velocimetry, so-called Lorentz force fiowmeters, it is 
necessary to predict the force and torque acting upon a complex- shaped magnet system 
by a turbulent liquid metal fiowing in pipes, ducts and open channels. Such predictions 
are usually performed by numerically solving the full set of three-dimensional equations 
of magnetohydrodynamics [3]. In order to be able to assess the reliability of such 
simulations it is necessary to have simple models that are amenable to rigorous analytic 
treatment. The goal of the present paper is to formulate and solve such a model. 

The model to be studied in the present work is a generalization of a previously 
studied problem [9] to the case when the orientation of a magnetic dipole is arbitrary. 
More precisely, the authors of |9] investigated the interaction between a moving plate 
and a magnetic dipole whose orientation is perpendicular to the surface of the plate. 
Here we relax the assumption about the orientation of the dipole and allow the dipole 
to be arbitrarily oriented. As will be shown in Section [21 this problem can still be solved 
exactly and all components of the torque and the force can be expressed analytically 
along with the electric and magnetic field. We illustrate our results in Section [3] using 
some representative plots of the three-dimensional structure of the eddy currents and 
the induced magnetic fields. In Section |4] we summarize our conclusions and indicate a 
possible illustrative application of the theory that would be interesting to investigate. 

2. The problem and its exact solution 

Consider a single magnetic dipole with dipole moment m located at a distance h above 
a fiuid or solid layer with electrical conductivity a and thickness w as shown in Fig. [1] 
The dipole is assumed to be arbitrary oriented in space. For brevity, if it is not specified 
otherwise, we assume that the index n denotes one of the Cartesian axes, m = me„. 
The general case can be considered as a linear combination of the three generic cases, 
n = x,y,z, any general orientation a can be expressed as a = m ■ {a^"^^}^, where 
{a^""^} = (a^^\ a^'^\ a^^^) is the tensor containing components for each n = x,y,z dipole 
orientation. The layer moves with an uniform velocity u = ue^ and extends in the z 
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Figure 1. Sketch of the problem: A plate with thickness w moves with constant 
velocity u and interacts with an arbitrarily oriented magnetic dipole. 

direction from to Zu- The thickness of the layer is then given hy w = \zu — Zd\- Here, 
the indices "ci " and "u " stand for "downward side" and "upward side", respectively. 
The borders of the layer in the x and y direction are supposed to be at infinity. This 
corresponds to the assumption that the distance h is much less than the horizontal 
size of the moving solid layer. We further assume that plate or the fluid moves slowly 
which is equivalent to the assumption that magnetic Reynolds number Rm = fiQauw 
is small. In magnetohydrodynamics (see e.g. [7j, [3]) this case is referred to as the 
quasistatic approximation and implies that the magnetic field associated with the eddy 
currents is much smaller than the applied magnetic field. The assumption Rm <^ 1 
makes the present problem amenable to rigorous analytic treatment. For instance, the 
magnetic Reynolds number of liquid aluminium (cr ^ 3 x 10^) flowing with a velocity 
u = 1 m/s in a pipe with a diameter of 0.1 m is approximately 0.3. In what follows 
we will use mainly Cartesian coordinates and will not apply an integration of Bessel 
functions in cylindrical coordinates as we did before in [9]. Our goal is to derive explicit 
compact algebraic expressions, which are useful for fast analytic calculations, for the 
primary magnetic field B(r), the electric potential 0(r), the eddy currents j(r) and the 
secondary magnetic field b(r). Then, all other quantities such as forces and torques 
acting upon the dipole can be readily written down. 

2.1. Primary magnetic field 

The first step in the solution of our problem is the specification of the magnetic field of 
the dipole which is referred to as the primary magnetic field. Let the magnetic dipole be 
at the point tq = (xq, yo, Zq) and we wish to calculate the magnetic field B(r) at a point 
r = {x,y,z). Before presenting formulae for the primary magnetic field, we introduce 
the following short notations of the vector R = r — tq: 

Ax = x-Xq, Ay = y-yo, Az = z - Zq, 
Ar^ = Ax^ + Ay^, R^ = Ar^ + Az^, 
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Now, the magnetic field of a magnetic dipole can be represented as [5]: 
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Tliese expressions appear more complicated than the formula normally used. However, 
the present notation will turn out to be convenient for the further analysis and it is the 
key to a compact representation of all fields. Then it follows, by taking ([T]) and m = me„ 
into account, that the components of the primary magnetic field can be conveniently 
given as: 
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where we introduced the abbreviation fi = nifio/Air, and is the partial derivative 
with respect to the coordinate specified by the direction n of the dipole. If the dipole m 
is arbitrary oriented in space, m = ma-e^ + myey + m^e^, then 9^ = m- V. For instance, 
for the dipole oriented parallel to the 2-axis, n = and n = z, and the components of 
the primary magnetic field are: 
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2.2. Electric potential 

The derivation of the electric potential is an extension of the method developed earlier 
in [9], and it is presented also in [6]. 

The electric potential (f) and eddy currents j are governed by Ohm's law for an 
electrically conducting material moving with velocity u, which can be written as: 

j = a(-V0 + uxB), V-j = 0. (5) 

(The second expression is not Ohm's law and does imposes the requirement that there is 
no source or sink of electric currents in the moving plate.) In general, the magnetic field 
in ([5]) is the sum of the primary magnetic field and the magnetic field associated with 
the eddy currents, referred to as the secondary magnetic field. As explained earlier, we 
assume that the magnetic Reynolds number Rm = fioavh, which characterizes the ratio 
between the secondary and primary field, is small, so we can approximate the magnetic 
field in Ohm's law by the primary field given by ([3]). 

The primary magnetic field is solenoidal, V x B = 0, and the velocity u is constant. 
By applying the divergence operator to Ohm's law (j5]) we therefore obtain the following 
equation for the electric potential: 

A</) = V • [u X B] = B ■ [V X u] - u ■ [V X B] = 0. (6) 
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Hence, is a harmonic function obeying Laplace's equation, A0 = 0. Any harmonic 
function in three-dimensional space can be expressed through derivatives or integrals of 
the function 1/R defined in ([T]), where the specific expression for depends on boundary 
conditions. In our case, the domain is infinitively large in x and y directions, and finite 
in the z direction. At the upper {z = Zy) and lower {z = Zd) surfaces of the plate. 



the vertical electric current must be zero, which is expressed as 



0, where the index s stands for either s = u oi s 



can directly express (9^0 by setting jz 
dz(p = uBy = ufi dydrA 



everywhere in the layer: 
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where the last equality directly follows from Now, if we change the order of 

differentiation and put dz as the first derivative, then the integration is trivial, and 
we come immediately to the final result: 

Az' 

where the second term (j)oo{x, y) appears as constant of the integration and so it does not 
depend on z. This term is determined from the additional boundary condition specifying 
that the electric potential vanishes when the dipole is at infinite distance from the layer. 
The result is: 
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where sgn(2;) = z/\z\ is the sign function with the property that sgn(A2;) = 1 for 2; > 2:0 
and sgn(A2;) = —1 for z < zq. One can readily verify that given by (|H]) is the 
solution of ([6]). To do this, we notice that = 0(l/i?), where 1/R is given by ([1]) and 
A[l/R] = A[tsinh~\Az/R)] = 0, hence, A0 = too. 

In Appendix A| it is shown that the components of electric current can be expressed 
through a stream function ip defined as j = V x (ipez). This stream function is given 
by: 
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In the short notation both and ip can be conveniently represented by means of the 
differentiation of the same auxiliary function $ as: 

0(") = V'^"^ = -u^crdA^, (11) 



$ = sgn(A2;) I tanh 



|A^| 



Inr 
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where an upper index n specifies explicitly the direction of the magnetic dipole. 

Because the electric currents are two-dimensional, the solution f|TT]) describes an 
electric potential and current stream function for any horizontal plane of the moving 



layer. Explicit formulae for 'ip and j are given in Appendix A These formulae can be 



immediately used for calculation without performing cumbersome differentiation. 

2.3. Secondary magnetic field and scalar potential 

The eddy currents j of the moving plate induce a magnetic field b which we refer to as 
the secondary magnetic field. When j is time-independent as in our case, the secondary 
magnetic field b(x, y, z) at a point r = (x, y, z) is uniquely determined by Ampere's law 

'^{x, y, z) = —V X b(a;, y, z). (12) 

On the other hand, j can be determined with the aid of the stream function which has 
been already above calculated. Eg. ( ITT]) : 

y, z) = Q{z)W X [e^^(x, y, z)] . 

Here the auxiliary step function 6(^) is introduced in order to ensure that eddy currents 
outside of the plate are zero. For instance, in the case of the layer of finite thickness, 
the function Q{z) is defined as: 



e(.) 



1 for Zd< z < Zu] 

for 2; < Zd and Zu < z, 



where Zd and Zu are the lower and upper borders of the plate. In the general case, the 
secondary field can be represented completely as a sum of the stream function and a 
gradient of a scalar function a in the form: 

b(x, y, z) = e{z) /io i^ix, y, z) + /ioVa(x, y, z), (13) 

Function 6(2;) enforces the condition that outside the plate the secondary field b is 
solely a gradient of the scalar potential a. This automatically ensures that V x b = 0, 
and then the problem how to calculate b is reduced into the problem how to calculate 
the scalar potential a. In order to accomplish this task one has to employ the constraint 
that b is divergence- free. By taking the divergence of (fT3|) . the divergence-free property, 
V ■ b = 0, results in the following equation: 

"1" 

'r 

where the second equality directly follows from Eq. ( ITT]) . So one has to solve now 
equation (IT^ with the boundary conditions | Va| — which describes that the secondary 
magnetic field vanishes at infinity. 

It is possible to solve Eq. f fTT|) in a straightforward way for any specific n = x,y, z. 
This requires explicit partial derivatives with respect to x and n and then cumbersome 
symbolic calculations for each separate index n by recruiting series expansion with Bessel 
functions and further integration. Such procedure has been performed in [S] for the 



Aa = —6(2;) dzi> = Q{z)ufiadxd„ 



(14) 
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particular case n = z, see also [Sj- However, there exists another technique to provide 
a general formula for arbitrary general n. For this we notice that a derivative of any 
function f{R), which depends on R solely, with respect to n = x,y,z can be taken as 
a derivative, with the opposite sign, with respect to Uq = Xo,yo,Zo, correspondingly, 
dnf{R) = —dnof{R), where R is defined in ([T]). In order to use this property of R we 
represent the scalar potential a by means of a generating (or primitive) function A twice 
differentiated with respect to the coordinates of the dipole, a = ufia dx^dnoA. We insert 
the above formula for a into (fT4|) and obtain the simpler Poisson's equation independent 
of n to determine A as: 

AA = e{z)^. (15) 

This equation can be readily solved without cumbersome integration and series 
expansion as we shall demonstrate below. Once we obtain the primitive function A, 
the secondary magnetic field can be represented by means of triple differentiation in the 
form: 

b(") = /io/i^wV [dx,dn,A] . (16) 

2.4- Scalar potential and secondary magnetic field for an infinitely thin layer 

Because the problem under consideration is linear, we treat the case of an infinitely thin 
layer separately. Once a solution for this case is found, all other cases can be presented 
through a superposition of thin layers. 

Let i be a location inside a thin moving layer with electric currents induced by 
an external magnetic dipole. Let us assume that the thickness of the layer is small 
compared to the distance between the dipole and the layer which is expressed by the 
condition w = \zu — z^l ^h. Finally let us introduce z = \zu + Zct\/2. In this case, the 
function Q{z), which is needed to set the position of the layer, must be taken in the 
limit Zd — )■ Zu- It gives Q{z) = w6{z — z), where w is the thickness of the layer which 
is inserted for the proper physical dimensionality. (Function Q{z) is dimensional in the 
case of the infinitely thin layer, but w will be omitted for other cases because there the 
proper dimensionality follows directly from an integration over the plate thickness.) 

Let Aq be the primitive function for the thin layer. Here the index is inserted 
to distinguish the present case from other cases. As follows from Eq. (fT5|) . everywhere 
in the space AAq = except for z = z. This means that Aq is a continuous harmonic 
function showing a discontinuity in the first derivative with respect to z ai z = z. As 



shown in Appendix B, the solution of f|T5|) is: 



Ao{z) = ^ { tanh"^ +\nr}, (17) 




R{zf =r' + \Az{z)\'\ \Az{z)\ = \z - zo\ + \z 

Notation |A2;(i)| looks ugly but this is because we want to stress an analogy in between 
functions R and R as explained below. In the given case of a infinitely thin layer, it 
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is superfluous to put argument z in the round brackets of the function Az, but this 
argument will be necessary afterwards, when we shall extend these results to the case 
of a layer with flnite thickness. 

Note, that Aq in (fT7|) and $ in (ITT!) are almost identical, that is, (w/2)$ corresponds 
to sgn(Az)y4o, with the mapping i? <H- and Az ■H- Az. Moreover, further analysis 
shows that this correspondence is valid not only for the entire functions Aq and but 
also for all their derivatives with respect to n = x,y,z and Hq = XQ,yo,Zo, that is, 
(w/2)9„$ O sgn{Az)dnAo, c?„$ = —9^0$, and dnAo = —dnoAo. The only exception 
is dzAo = dzgAo because {z — Zq) is taken in Aq as an absolute value. By taking 
this property into account and recalling that both the eddy currents j and secondary 
magnetic fleld b can be represented through a triple differentiation as: 

^ -V X [d^d^^ e„ — ^ = V [d,,dn<Ao] ■ 



Therefore without computing derivatives, it is possible to build up the following 
correspondence between the secondary magnetic fleld and eddy currents: 

sgnjA^) .^^-^ ^ ^ sgn(A2;) 

2 

sg n(A2:) .(^^ , ) sgn(Az) 
~2 ^ y ' y 2 
(The difference in sign for n = z is because 9^$ = —d^^,^ while O^Aq = OzqAq.) The 
remaining ^-components of the secondary fleld are: 



^ forn = x,, (18) 
bi:) ^ forn = .. (19) 



^,)^sgn^^.(, 
^(,)^sgn^^^,,)_^,,)^_ (20) 

The last relationship for fei^^ is the conjecture of the fact that the $ is a harmonic 
function, i.e. A$ = 0, hence, {jy^^ — ji^^)/{j2au) = {d^d^ + dydy) ^ = [—d^dzdz^] 
and {w/2) d^d^dz^ O sgn(A2;) d^dzdzAQ. According to the mapping fll8tl20p . the tensor 
of the secondary magnetic fleld can be now written down by means of the tensor 
{j^")} given in ( lA.lOl) for eddy currents. This tensor is presented in Appendix C and 
can be straightforwardly used in algebraic calculations in contrast to the implicit results 
published by Priede [6]. 



2. 5. Secondary magnetic field for half-space and for the layer with finite thickness 

Because the problem under consideration is linear, the secondary fleld of any moving 
plate can be represented by summing all the flelds originating from the layers of the 
plate derived in the preceding section. Mathematically, this means an integration of the 
secondary magnetic fleld tensor, cf. Eq. ( ]C.1|) . over the thickness of the plate. The limits 
of the integration in the case of half-space are z^ and z^ — )■ — oo. Then, the integration 
should be carried out for the Pk coefficients only because the functions fk and gk do not 
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depend on z. By performing this procedure formally, we obtain 

Poo,k = I Po,k{^) dz, (21) 

J — oo 

where Po,fc and Poo,fc are tensor coefficients of an infinitely thin layer and half-space, 



correspondingly. Explicit formulae for Poo.fc are presented in Appendix D 

When the plate is of finite thickness, < z < Zu, the coefficients Pk are obviously 
expressed through the coefficients for half-space in the form: 

Po,k{z) dz= Po,k{^) dz- Po,k{^) dz 

(J J —oo J —oo 

= Poo,k{Zu) — Poo,k{Zd), (22) 

where Poo,k{zu) and Poo,k{zd) are computed according to f ID.ip with Az and R taken at 
points Zu and z^. 

2. 6. Magnetic potential energy, force and torque 

In the general case, the potential energy [/ of a magnetic dipole m in a magnetic field 
B isU = —m ■ B, see e.g. In our case, the magnetic field is the induced secondary 
magnetic field b = nofxau [m ■ {6*^"^}-^ ■ m] , hence, the potential energy of the dipole is 
written down as: 

U = - fionau [m ■ {6^"^}^ ■ m] . 

= - /io/icTM [ml bi^^ + ml bf + ml bf + 2 myb^y^\ , (23) 

where the last relationship is obtained by taking into account a symmetry of the 
secondary field tensor b^y^ = b^x \ by^^ = —b^z \ and b't^ = —b^^\ 

Force -F acting upon the dipole can be calculated as a negative gradient of the 
potential energy: 

F= - VU, F, = fiofiau [m • {f^Y " mj , (24) 

where a force field tensor {f is a partial derivative of the secondary magnetic field 
tensor with respect to the coordinate axis i = x,y, z, i.e.: 

{/i")} = 9.{6(")}. (25) 

Because the tensor {b^"'^} has been found above, algebraic formulae for the tensor 

can be computed as well. They are presented by Eq. ( ]E.mE.3p . Notice, that these tensors 



contain 27 terms (3 times 9), while due to symmetry and algebraic simplifications they 
can be given by means of linear combinations of five coefficients Qk, k = 0, . . . ,4. 
The torque T acting upon the dipole is: 

T = mxb = -fiofiau (m x [{b^")}^ ■ m]) . (26) 

It can be computed also also by means of the tensor {b*^")} presented explicitly in (IC.ip . 

Above results for the potential energy, force and torque are of the same structure 
when the secondary magnetic field is induced by a discrete ensemble of magnetic dipoles 
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or by a continuous macroscopic magnet. (The latter can be obtained by an integration 
over the space occupied by the magnet.) However in our simplest case, in order to 
calculate the force and torque acting upon a single dipole, we just put r = and z = zq 
in all the tensors found already. Mathematically, these tensors are sums of the terms, 
which have functions fk and gk as factors with coefficients Pk or Qk- The functions fk 
and gk are equivalent to cos kcj) and sin kcj), correspondingly, here, tan0 = Ay/ Ax, and, 
hence, fk and gk are zero at r = for k > 1. Hence, at r = 0, all the tensor terms 
having factors Pk or Qk with k > 1 must vanish due to their products with fk and gk- 
This greatly simplifies all the formulae, and we obtain finally: 

b = fio^au Pol.^Q .^^^^ {-m^e^ + m^e^} , (27) 

F = fiofiau Qo I, .=0,2=^0 + + 4m^)ea; + 2m^myey} , (28) 

T = fiQfiau Polr.=o,2=2o {i^xfriye^ - {m\ + m^)ey + m^^m^ej , (29) 

where the non-zero tensor coefficients Pq and Qo at r = and z = zq are presented in 
Appendix F in formulae (IF.lt IF.2P for a thin layer, (IF.3t IF.4P for a half space, and ( IF.5[ 



IF.6P for the layer with finite thickness. 

3. Spatial structure of eddy currents and secondary magnetic field 

In order to illustrate the analytical results of previous sections, we plotted eddy currents 
and the secondary magnetic field by using the formulae derived above, Eq. IA.6P and 



( ]C.1|) . The eddy currents below are shown as contour lines of the eddy-current stream 
function given by Eq. (lA.6p . Solid lines correspond to positive and indicate 



counterclockwise direction of the eddy current; whereas dotted lines are for negative 
contour levels and imply clockwise direction of the eddy current. The secondary 
magnetic field is given in the form of selected three-dimensional magnetic field lines 
computed by using Eq. ( IC.ll) . According to the right-hand rule, the secondary magnetic 
field lines emanate from the solid positive loops and penetrate the negative dotted loops. 

Figure |2] shows behavior at niy = and continuously changing nix and rUz provided 
that -|- niy + = 1. This corresponds to the case when the dipole is confined 
to the plane spanned by the velocity vector and the direction is normal to the surface 



of the plate. The results of Fig. 2(a) have already been given in [9]. In this case, 
= my = 0,mz = 1, and there are two loops of eddy currents whose centers are 
arranged parallel to the direction of the plate movement. The secondary magnetic field 
lines form a cage. When contribution increases, both loops deform and the positive 



loop grows in its size and shifts under the dipole as shown in Fig. 2(b) Then, the case 



of Fig. 2(c) corresponds to a magnetic dipole, which is oriented along the x-direction, 
i.e. = 1, my = 0, m^ = 0, where the positive eddy current loop is located completely 
under the dipole, and two negative loops are at the periphery. As a result, the three- 
dimensional lines of the secondary magnetic field are either completely in the y — z 
plane, if they start at x = 0, or the lines are closed at the peripheral negative loops if 



they start at y = 0. This behavior is emphasized by Fig. 2(d) which refers to the same 
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parameters as Fig. 2(c) but plotted from a different perspective. Then, Fig. 2{e] 



IS a 



mirror of Fig. 2(b) with the difference that the directions of the eddy currents and of 
the secondary magnetic field are opposite. The same is true for Fig. 2(f) as compared 



to Fig. 2 (a; 



Figure [3] shows the behavior at = and continuously changing rUx and 
provided that + + = 1. This corresponds to the case when the magnetic 



dipole is located in a plane parallel to the surface of the layer. The case of Fig. 3 (a 
0, niz = is the same as Fig. 2(c) and is repeated here in order to guide 



m 



1, m,, 



the eye. By rotating the dipole, the central loop turns together with the lines of the 
secondary magnetic field. Fig. 3(b) Then, when the main axis of the central loop is 
parallel to the diagonal of the plate, it splits into two equal positive loops, and thus one 
observes four alternating eddy currents loops for = 0, niy = 1 = 0, see Fig. 3(c 



The magnetic lines emanate from the positive loop and go into the neighboring negative 



loop. If one looks at the system from above, as in Fig. 3(d) being the same as Fig. 3(c) 
but differently oriented, one can see that the line r = 0, where the dipole is located, is 
the place where the magnetic field lines meet and turn out, holding the secondary field 
equal always zero. Then, the Fig. 3(e) is a mirror of Fig. 3(b) and Fig. 3(f) is a mirror 



of Fig. 3(a) with the difference that directions of eddy currents and secondary field are 
opposite in respect to each other. 



4. Conclusion 



An analytic solution is obtained for the electromagnetic interaction between an 
arbitrarily oriented magnetic dipole and a moving electrically conducting solid layer. 
The solution includes the electric potential and eddy electric currents inside the layer 
as well as the scalar potential for the secondary magnetic field outside the layer. The 
formulae obtained are useful for a Lorentz force velocimetry, because a moving solid 
plate can be considered as the limiting case of the mean velocity profile of a turbulent 
liquid metal flow when the Reynolds number (Re) tends to infinity. It has been shown in 
[9] that the signal of a Lorentz force velocimetry system interacting with a turbulent pipe 
flow converges towards the signal generated by a moving solid cylinder when Re— > oo. 
Hence the present model provides a benchmark for numerical simulations of Lorentz 
force velocimetry in channel flows at high Reynolds numbers. 

It should be finally mentioned that the present model could be applied in a 
straightforward way to describe the behavior of an educational experiment which 
consists of a spherical permanent magnet rolling down an inclined aluminium plate. 
Depending on the initial orientation of the magnetization the sphere undergoes different 
kinds of motion including straight motion with periodically changing velocity, wiggling 
motion and irregular motion. With minor changes the present model could be extended 
to cover this case. 
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(m^ rriy m^) = (0.00 0.00 1.00) 



(m m m ) = (0.71 0.00 0.71) 

^ X y z' * ' 




(a) 



(b) 



(m m m ) = (1.00 0.00 0.00) 

^ X y z' * ' 




(c) 



(d) 



(m m m ) = (0.71 0.00-0.71) 

* X y z' * ' 



(m^mymj = (0.00 0.00 -1.00) 




(e) 



(f) 



Figure 2. Eddy currents and secondary magnetic field for the case m.y ~ 0, i.e. when 
the magnetic dipole is located in the plane spanned by the velocity vector and the 
direction normal to the surface of the plate. The points where the secondary magnetic 
field lines emanate are selected manually inside and nearby the eddy-current loops in 
order to guide the eye. 



Interaction of a magnetic dipole with a slowly moving electrically conducting plate 14 



(m^ % m^) = (1 -00 0.00 0.00) m m^) = (0.71 0.71 0.00) 




(c) (f) 



Figure 3. Eddy currents and secondary magnetic field for the case = 0, i.e. when 
the magnetic dipole is located in a plane parallel to the surface of the moving plate. 
The points where the secondary magnetic field lines emanate are selected manually 
inside and nearby the eddy-current loops in order to guide the eye. 
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Appendix A. Eddy currents expressed through a stream function 

By inserting (j){x, y, z) from ([8]) into ([5]), we find tlie components of electric current as: 



, -1 

tann — 
R 



(A.l) 



iy = ^ {-dy<P - u B^} = ufia ^d^d^dn tanli ^ ^ - dyCpoo^ , (A.2) 
3z = 0, (A.3) 
To simplify the above expression for j^, we used from (|3]) and the fact that 



1 -1 

tanh 

R 



idyy + 



, -1 Az' 
tanh — 
R 



(A.4) 



The obtained eddy currents are purely horizontal. Due to their two-dimensionality they 
can be expressed as 

j = V X (tpe,) 

where the stream function ip is 



ilj{x, y,z) = - una ^yd^dn 
ipoo{x,y) = - lim dccdn 



, -1 ^z' 
tanh —— 
R 

tanh"'^ — - 
R 



= sgn( Az) dxdn In r. 



(A.5) 



The function ipooix, y) is a conjecture of the function (pooi^, y) and it is needed to nullify 
eddy currents at the infinite distance from the dipole. 

When the dipole is oriented arbitrarily, the general expression for ip can be given 
explicitly with the aid of the tensor {■i/?^"')} = as follows: 

= u/ia (m ■ {V'^^^D , {V'^"^}= ( i^o + i^2/2 F2g2 F,h) , (A6) 
Ax 



where 



9i 



r 

Ay 



fk — /i/fc-1 — gigk-i, 
9k = gifk-i + /i/fc-i; 



and 



Az 



Az sgn(A2;) 



|A^| 



(A.7) 
(A.8) 

(A.9) 



The eddy currents j = ujia {m • {j*-"^}"^) can be represented with the aid of the tensor 
-which is: 

/ Ci^l + 6*3^3 C1/1-C3/3 0292 \ 

-SCh-C^fs -Cigi-C,gs Co - C2/2 , (A.IO) 
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where 



Co = 777^ -7^, Ci = sgn(Az)^-^, = (A.ll) 



Cs = sgn(A2) 



3r 1' 
IR2 + ~ 



2 



lA^I 



- 1 



(A.12) 



R 

The coefficients Ck tend to zero when | A2;| goes to infinity. To find the eddy currents 
under the dipole, r = 0, one has to make a series expansion around r = 0. It gives: 

CoUo-^j-l+IAT^I+OM, (A.13) 

o 2 r f- 2 ^ 

C,U„^{-5I1^}^0M (A.16) 

Thus at r = 0, Co is not zero. It is negative, hence, under the dipole, the direction of the 
eddy current is opposite to the y-axis, when the dipole is oriented along the 2;-direction 
in agreement with the right-hand rule. These series expansions results will also be useful 
afterwards in the analysis of the secondary magnetic field and forces acting upon the 
dipole. 

It is important to notice that the functions fk and gk above depend on x and 
y only and do not depend on z. This will be be useful below when we compute 
secondary magnetic field by integrating over the plate thickness Zd < z < Zu ■ Moreover, 
these functions fk and gk are specially determined is such a way, that if the origin of 
the Cartesian coordinates is selected on the line passing through the magnetic dipole 
perpendicular to the plate, xq = yo = 0, then in the cylindrical coordinates, tan ip = y/x, 
fk = cos kip and gk = sin kip, i.e. the formulae above represent Fourier decomposition. 
There must be no angular dependence in a series expansion around r = precisely at 
the point r = 0, therefore, Cq is the only non vanishing coefficient in (lA.lSp . 

Appendix B. Primitive function for a thin layer 



Here, we prove that the solution of f[T5|) for an infinitely thin layer is given by: 

Az{i 
R{z 

where 



Ao{z) = - { tanh-i + in r J> , (B.l) 



R(z)'^ = + |Az(z)|^ \Az{z)\ = \z - zo\ + \z - z\. 

Function |A2;(i)| can be understood from the following geometric interpretation. It is 
the sum of two terms, the first term, |i — ^ol? is a distance between zq and z, and the 
second term \z — 2;| is a distance between z and z. Hence, the layer located at z = z 
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\Az(z) 



subdivides the whole space into two half-spaces, z > z and z < z, therefore the function 
Az{z) in straight brackets takes the following form: 

\z + Zo-2z\ a {z - z){z - zo) > 0; 
\z - zo\ = \Az\ if {z - z){z - zo) < 0. 

The first line is applicable when z and zq are both in the same half-space, hence, 
Az{z) = z + zq — 2z in this case; the second line is applicable when z and zq are 
in different half-spaces, hence, Az{z) = z — Zq in this case. Then, it is important that 
|A2;(5)| shows a discontinuity in the first derivative ai z = z, i. e.: 

+ 1 if 2; > z, 
— 1 if 2 < z. 



1 d^ 

hence - \Az{z)\ = 6{z 



Now it is easy to ascertain by means of straightforward calculations that AAq = 
everywhere except for the point a.t z = z while at this point the formula (fT5|) is valid 
because: 




w \Az{ 


5)1 


2 i?3 ' 


\Az{z)\ 


'd 




dz 


d'Ao^ 


d^A 




dz"^ 



+ 



\Az{z 

w 6{z — z 



1 rf2 



\Az(z) 



Rdz^ 
1 

R' 

where the equality R = R a.t z = z is taken into account. 

Appendix C. Tensor for the secondary field of a thin layer 

As follows from the mapping fll8tl20p between induced field and eddy currents, and the 
tensor (1A.10|) for the eddy currents, the tensor of the induced field can be written down 
as: 

/ -3/1P1-/3P3 -giPi-gsP^ /2P2-P0 



(C.l) 



-(7iPi - gsPs fsPs - hPi g2P2 

\ -^72^2 -4/1 Pi 

where the coefficients Pk = w Po^k for an infinitely thin layer are obtained from the 
mapping Po k ^ sgn{Az)Ck/2, /c = 1, . . . , 3, by keeping in mind that R and Az are to 
be used instead of R and Az. Specifically, 



-Po,o 
Po,2 



sgn(Az 
—sgn{Az 



3r' 



4P2 2P3 
3r2 



P 



^,1 
\Az 



0,3 



3r| A2I 
8P5 

3r r 
r 



(C.2) 



4P2 



1 



\Az\ 



R 



4R5 ' 2 r3 

The series expansion flA.13p -f lA.16p is also applicable for f lC.2p . One can see that at 
r = Pq.o is not vanishing only, therefore, as follows from f lC.ip . the nonzero tensor 



components at r = are bx 



only. 



Interaction of a magnetic dipole with a slowly moving electrically conducting plate 18 



Appendix D. Secondary magnetic field tensor for a half-space and for a 
layer with finite thickness 

For a half-space and for a layer with finite thickness, the secondary field tensor is of the 
same structure as for an infinitely thin layer, Eq. ( IC.ll) . with the following coefficients 
for the half space: 



00,0 



00,1 



Az 

r 

16^3' 



00,2 



00,3 



sgn(A2; 




\Az\ 



R 



(D.i; 



1 sgn(A^)i? 



IQR^ ARr 



2r3 



R 



where Az and R should be computed at Zu- 

For the layer with finite thickness, the tensor coefficients are: 

Pk{Zu-, Zd) = Poo,k{^u) — Poo,k{^d), k = . . . A] 



(D.2) 



where Poo,ki^u) and PooA^d) are computed according to fID.ip with Az and R taken at 
points Zu and Zd- 



v 

/ 

\ 
/ 

v 



Appendix E. Force field tensor 

The force field tensor, {/f"^} = di{b'^''^ can be written down in the following way: 

f 3Qo + 2f2Q2 + /4Q4 92Q2 + 94Q4 'ifiQi + /3Q3 \ 

g2Q2 + giQi Qo - UQi giQi + g-sQs (E.i) 
-3/1Q1-/3Q3 -giQi-g3Q3 ^Qo + '2f2Q2 J 

fl'2<52 + 5'4Q4 Qo - f4.Q4 giQi + g-sQs ^ 
Qo - /4Q4 fi'2<52 - giQi fiQi - fsQs , (E.2) 
-giQi - gzQ-i hQz - hQi '2g2Q2 J 

-SfiQi - hQ^ -giQi - g^Q^ AQ^ + 2f2Q2 \ 
-giQi-g^Q^ hQ^- hQi 2g2Q2 ■ (E.3) 

-AQo-2f2Q2 -2(72^2 -4/iQi / 

The coefficients Qk = w Qo^k for a infinitely thin layer are: 



Qo,o 
Qo,2 
QoA 



3\Az\ ( 5r 



8R^ 1 2i?2 
ISr^lA^I 
8^7 ' 



Qo,i 

Qo,3 = 



= sgn(A2; 
sgn(A2 



3r f 5r 



2i?5 14R2 



8i?7 



(E.4) 



3|Az| 


" 5^2 


3 1 


3 


lA^I 






h ^ + ^ 




' . ' 1 


2^3 


.8^4 


4/22 r2 




R 



For a half-space, the coefficients Qk = Qoo,k are obtained by using the same method as 
done before for the coefficients Poo,fc- This involves an integration over the thickness of 
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the plate, Qoo,k = /f^ Qo,k{z) dz, and the results of the integration are: 



Qoofl 



16i?3 [2i?2 



3r2 



Qoo,l 


— sgn(A2;^ 


3r A2; 
16^5 ' 


Qoo,2 


3r2 
16^5' 




Qoo,Z 


— sgn(Az' 


1 i^fi 


1 + ^ 
4i?2 


1 

+ ^ 





00,4 



2 4 



+ 



3i? 
2^ 



- 1 



(E.5) 
(E.6) 
(E.7) 

(E.8) 



For a plate with finite thickness we have: 

Qk = Qoo, k{zu) - Q 00, k{zd), A; = 0...4. 

Appendix F. Non-zero tensor coefficients at r = and z = zq 

The coefficients for a thin layer of thickness w <^ \zq — z\ located at z are: 
PI -PI _ sgn(2;o - z) 



Qo\r=0,: 



wQo,ol=o,z- 



lQ\zo — z\ 
3 



-w- 



(F.l) 
(F.2) 



The same coefficients for a half-space are obtained by integration over z from Zu up to 



infinity 



PI — P \ 

^ 0\r=0,z=zo — ^ oo,0|,,=o,z=2o 



1 



Qoo,o| 



32 \zq — 

1 



,,z=zo ^^,ulr=0,z=zo 128\zo-Zu\^' 

Finally, for a plate with finite thickness, Zd < z < Zu'- 



Pr 



Olr=0,2= 



1 



Qo\r=0. 



^=^° 128 \\zo-Zd\^ \zo-Zu\-' 



(F.3) 
(F.4) 

(F.5) 
(F.6) 
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